perm filename MODD.F4[JC,MUS] blob
sn#007329 filedate 1972-07-16 generic text, type T, neo UTF8
00100 DIMENSION W1(512),FUNC(512),SU(50),ST(500)
00200 GO TO 60
00300 40 CALL HYDPOG(2)
00400 60 TYPE 50
00500 50 FORMAT(' TYPE A,CF,MF,I1,I2,STEPS'/)
00600 ACCEPT 100,A,FC,FM,XI1,XI2,STEPS
00700 100 FORMAT (6F)
00800 CALL TYPLOC(-300,-512)
00900 CALL DPYSET(1,SU,50)
01000 CALL DPYBRT(1)
01100 CALL AIVECT(0,0)
01200 CALL ALINE(-264,0,256,0)
01300 CALL ALINE(-256,-256,-256,256)
01400 CALL DPYOUT(1)
01600 XI=XI1
01700 XI1=(XI2-XI1)/STEPS
01800 N=STEPS+1
01900 DO 400 N1=1,N
02150 I=-1
02175 CALL DPYSET(2,ST,500)
02187 CALL DPYBRT(5)
02200 DO 200 J=1,512
02300 X=J-1
02400 Y=6.2832/512.
02500 FUNC(J)=A*(SIN((FC*X)*Y-XI*(COS((FM*X)*Y))+XI))
02600 IY2=FUNC(J)*256.
02700 IF(I.EQ.0)GO TO 150
02800 CALL AIVECT(-256,IY2)
02900 I=0
03000 GO TO 160
03100 150 CALL SVECT(1,IY2-IY)
03200 160 IY=IY2
03300 200 CONTINUE
03400 CALL DPYOUT(2)
03500 XI=XI+XI1
03600 400 CONTINUE
03650 TYPE 403
03675 403 FORMAT(' TYPE 0 TO CONTINUE'/)
03700 ACCEPT 300,M
03800 300 FORMAT (I)
03900 IF(I.EQ.0)GO TO 40
04000 END